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Abstract 

Aims. We investigate the surface nitrogen content of the six magnetic O stars known to date as well as of the early B-type star r 
Sco. We compare these abundances to predictions of evolutionary models to isolate the effects of magnetic field on the transport of 
elements in stellar interiors. 

Methods. We conduct a quantitative spectroscopic analysis of the sample stars with state-of-the-art atmosphere models. We rely on 
high signal-to-noise ratio, high resolution optical spectra obtained with ESPADONS at CFHT and NARVAL at TBL. Atmosphere 
models and synthetic spectra are computed with the code CMFGEN. Values of N/H together with their uncertainties are determined 
and compared to predictions of evolutionary models. 

Results. We find that the magnetic stars can be divided into two groups: one with stars displaying no N enrichment (one object); and 
one with stars most likely showing extra N enrichment (5 objects). For one star (0' Ori C) no robust conclusion can be drawn due 
to its young age. The star with no N enrichment is the one with the weakest magnetic field, possibly of dynamo origin. It might be a 
star having experienced strong magnetic braking under the condition of solid body rotation, but its rotational velocity is still relatively 
large. The five stars with high N content were probably slow rotators on the zero age main sequence, but they have surface N/H 
typical of normal O stars, indicating that the presence of a (probably fossil) magnetic field leads to extra enrichment. These stars may 
have a strong differential rotation inducing shear mixing. Our results should be viewed as a basis on which new theoretical simulations 
can rely to better understand the effect of magnetism on the evolution of massive stars. 

Key words. Stars: massive - Stars: atmospheres - Stars: fundamental parameters - Stars: abundances - Stars: magnetic field 



1. Introduction 

The evolution of massive star s is governed by two m ain physi- 
cal processes: their mass loss (Chiosi & Maede r'1986) and their 
rotation (Maeder & Meynet 2000). All additional processes di- 
rectly affecting these properties are thus major contributors to 
' the fate of massive stars. This is the ca se of magnetism. Since the 
' pioneering work o f iBabel & Montme rle ( 1997), it is known that 
the presence of a large scale surface magnetic field will modify 
massive stars winds. By deflecting material along the field lines, 
" ma gnetism affects the wind geo m etry and structur e. Simulations 
! by lud-Doula & Owockil (l2002h : lUd-Doula et all (l2008l l2009h 
have shown that equatorial overdensities can be created. Material 
can either fall back onto the stellar surface or be ejected in the 
magnetic equatorial plane depending on the interplay with ro- 
tation. In the case of very strong magnetic fields, observations 
and simulations show that a rigid magnetosphere strongly affects 
the wind structure, creating caps of stel lar material in the min - 
ima of the magneto-rotational potential dTownsend et al.ll2007h . 
Since the wind flow is deflected by the presence of a magnetic 
field, the way the star loses its angular momentum is affected. In 
addition, the internal structure of the star is also affected, forc- 



* Based on observations collected at the CFHT and the Telescope 
Bernard Lyot. 



ing the star to follow new evolutionary paths. Given the rele- 
vance of massive star evolution to many fields of astrophysics 
(galactic chemical evolution, interstellar dynamics, supernovae 
and gamma-ray bursts, nucleosynthesis...), it is crucial to under- 
stand the magnetic properties of these objects in order to con- 
strain their influence on the evolution and fate of OB stars. 

The recent development of powerful spectropolarimeters 
has revolutionized our view of massive stars magnetism. The 
first detection o f a magnetic field on the Trapezium O-type 
star 0' Ori C (iDonatiet al.1 l2002h dates back to less than 
a decade. This spectroscopically varying and strong X-ray 
emitting star has revealed a surface kilogauss field organized 
mostly as a dipole. Subse quent de tections have b een made on 
HD 191 612 (.Donati et aP l2006al). H P 148937 (iHubrig et al.1 
2008; W ade et al.ll2011bl). H P 57682 rminhut et al.ll2009h and 
HP 108 (Martins et al.ll2 01(y). Field strengths are of several hun- 
dred Gauss to a few kilogauss and the geometry, when con- 
strained, is consisten t with a dipole. T he only exception is the 
O supergiant ^ Ori A dBouret et al.ll200 8) for which a more com- 
plex field topology was tentatively derived, together with a re- 
markably weak field of 50-100 G. This raises the question of 
the origin of magnetism in massive stars. The debate is ongoing. 
A fossil origin, stable over Myrs and wit h a rather simple ge- 
ometry, is supported by some simulations (Braithwaite & Spruit! 
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l2004t iDuez et alJl2010h . In that case, the star retains the orig- 
inal magnetic field of its parental molecular cloud or dynamo 
field of its convective progenitor. Alternatively, dynamo pro- 
cesses may be at work. In analogy with low mass stars where 
a convective envelope (or even a fully convective structure) ex- 
ists, magnetic field might be generated in the convective core 
of high mass stars and subsequently transported to the surface 
jciiar bonneau & MacGregor 2001; Brun et al. 2005). A dynamo 
may also operate in the radiative z ones of massive stars (Spruit 
l2002HMacDonald & Mullanl2004 l). The field would be produced 
by the shear in this region. This idea has be en tested in numer- 
ical simulations, with contradictory results jBraithwaitd 120061; 
IZahn et al.ll2007h . 

However, assuming that this dynamo process was able to 
produce and maintain a magnetic field, evolutionary calcula- 
tions have been conducted to test the effects of this type of 
mag netism on the evolution of massive stars ( Maeder & Mevnetl 
l2004l) . They find that the internal structure is strongly affected, 
with a completely different rotation profile in the star compared 
to non-magnetic calculations including rotation. A consequence 
of this re-organization of the stellar interior is a modification of 
the transport of angular momentum and of chemical elements. 
In particular, iMevnet & MaedeJ (l2005h studied the effects on 
the surface chemical appearance of the star. They show that the 
presence of a dynamo generated magnetic field significantly in- 
creases the efficiency of mixing, leading to strong nitrogen and 
even helium surface enhancement. Further analysis of the feed- 
back effects of magnetic braking balanced these results though, 
highlighting the importance of the type of rotation profile (dif- 
ferential versus solid-body) in the final chemical appearance of 
the star (Meynet et al. 2011). 

Although still limited, the sample of magnetic O stars is now 
populous enough to test these predictions. An important side ef- 
fect of the spectropolarimetric observations is the production of 
very high S/N ratio optical spectra at high resolution. These are 
perfectly suited for abundance studies. In this paper, we thus an- 
alyze the surface nitrogen content of the magnetic O stars known 
to date as well as of the B0.2V star r Sco. Our aim is to see if 
peculiar chemical patterns are observed and if so, to what extent 
they can be explained by the current generation of evolutionary 
models. The goal is to identify specific patterns that can only 
be attributed to the presence of a magnetic field. This is par- 
ticularly important in the context of explaining the population of 
N-rich slowly rotating stars in the Magellanic Cloud reported by 
iHunter et al.i (I2008ll2009l) . Such stars are not explained by rotat- 
ing evolutionary tracks and have been suggested to be magnetic. 
The observations are presented in Sect.|2l We describe the analy- 
sis with atmosphere models in Sect. [3] The results are discussed 
in Sect. ID We show how the chemical properties of the magnetic 
O stars can be interpreted in terms of evolutionary models with 
rotation and with or without magnetic field. We finally summa- 
rize our conclusions in Sect. |5] 

2. Stellar sample and observational data 

2.1. Observations 

The observations were collected with echelle spectropolarime- 
ters ESPaDOnS at CFHT (Canada-France-Hawaii Telescope) 
and NARVAL at T BL (Telescope Ber nard Lyot), partly from 
the MiMeS project dWade et al.ll201 l a^. The spectra cover the 
wavelength range 3700-10500 A with a resolution of 68000 
(ESPaDOnS) and 65000 (NARVAL). We gathered multiple se- 
quences for each of the sample stars, each sequence consisting 



of four individual subexposures taken in different polarimeter 
configurations. From these four individual subexposures, mean 
Stokes parameters were derived. All frames were processed us- 
ing Libre-ESpRJT (iDonati et al.ll 19971) . a fully automatic reduc- 
tion package installed at CFHT and TBL for optimal extrac- 
tion of ESPaDOnS and NARVAL spectra. The peak signal-to- 
noise ratios (SNRs) per 2.6 km s"' velocity bin depends on the 
instrument, on the exposure time and on weather conditions. 
Complete references to the spectropolarimetric data are gathered 
in Table lA.fl A summary of the data available is given in Tab.[T] 
We used the average of the entire set of spectra for each star, en- 
suring that no variability in the key diagnostic lines was present 
among the data. 

2.2. Sample Stars 

Our sample is composed of the six magnetic O stars known to 
date. Here we give a few details about each of these objects. The 
journal of observations is shown in Table lA.fl 

HD 191612 is a member of the Of?p class, first intro- 
duced by Walb orn (1972). Such stars exhibit, among other pe- 
cuharities, variabilit y in their Balmer, He i, C iii and Si iii lines. 
iDonati et al.l (l2006al ) reported the detection of a kilogauss mag- 
netic field on its surface. It is presently well-established that 
HD 191612 displays two re current sta tes, with a p eriod of ap- 
proximatelv 538 davs (Walborn et all20 03: Walbor n & HowarthI 
120071; iHowarth et al II2007I) . These two states are referenced here 
as "quiet" (i.e. 08 spectral type. Ha in absorption) and "active" 
(i.e. 06 spectral type, strong wind signatures). We decided to 
focus on the quiet phase since the photospheric lines, used to 
derive the stellar parameters and surface abundances, are less 
contaminated by wind emission than in the active phase. We re- 
lied on data taken on August 23'^'' 2009, corresponding to a phase 
of 0.393 (see Wade et al. 2011c). The journal of observations in 
Table [llpresents additional information. 

HD 108 is the second Of?p star of our sample. It is very 
similar to HD 191612, with a long-term variability both photo- 
metrically and spectroscopically. However the timescale of this 
variability is much longer and is estimated to be of about 55 
years. The He i and Balmer lines are the most affected, leading 
to a spectral type changing from 04 to 08. |]VIartins et al.l ( l2010h 
reported the detection of a (most likely) kilogauss magnetic field 
on its surface. We used the average of the July- August 2009 sub- 
set of the Martins et al. observations in the present study since, 
as for HD 191612, this coiTesponds to the minimum of extra 
emission in the Balmer and He i lines. 

HD 148937 is the last member of the Of?p class in our sam- 
ple. Its variability is much more limited (lower amplitud e) than 
that ob served in HD 191612 and HD 108 spectra ( Naze et al.l 
l2008bh . The periodicity of the variations of about 7 days, is ex- 
tremely short compared to the long periods observed in the two 
other Of?p stars. The HD 148937 spectra were collected during 
six nights in June 2010. 

0' Ori C (HD 37022) is a mid-late type O d warf. It is the 
first O star for which a magnetic field was detected (IDonati et al.l 
I2OO2 I). It is known for its variability, with a period of approxi- 
mately 15.4 days (Stahl et al. 1996; Wade et al. 2006), the most 
striking variable spectral features being an inverted P-Cygni pro- 
file in the Heii/14686 line and strong emission features in the 
Balmer lines (especially Ho- and HjS). The equivalent width of 
some He and metal lines also vary as a function of the phase 
(Simon-Diaz et al. 2006, and references therein). 

Spectropolarimetric observations of f Ori A were collected 
with NARVAL in 2007 October. Altogether, 292 circular polar- 
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Table 1. Summary of observations. 



Star 


Tn^triimpnt 


Tof al Pvn tirnp 
i.yji.aL J_iA.u. Liiiic 


N^iirnhpr of cnppfrM 


tvnipnl npr *;npptriirn 


HD 108 


NARVAL 


34 h 20 min 


37 


150-500 


HD 148937 


ESPADONS 


4 h 00 min 


6 


200-700 


HD 191612 


ESPADONS 


2 h 40 min 


8 


100-350 


0' Ori C 


NARVAL 


53 min 


2 


500-1300 


^Ori A 


NARVAL 


13 h 52 min 


336 


1000-1400 


HD 57682 


ESPADONS 


3 h 4 min 


5 


400-800 


T ScO 


ESPADONS 


13 min 


5 


600-1200 



Notes. The complete journal of observations is given in Appendix lA.il 



ization sequences were taken. The peak SNR per 2.6 km s ' ve- 
locity bin varies between 800 and 1500. f Ori A probably has a 
short rotation period of approximately 7 days and a weak mag- 
netic field (Boure t et al. 2008). Although some spectral lines ap- 
pear variable, the main diagnostics used to derive photospheric 
parameters remain stable during the rotation cycle. 

HD 57682 is an 09 sub-giant. Its magnetic field of strength 
1.7 kG was reported by G runhut et alj (l2009b . HD 57682 dis- 
plays spectral variability but on a much lower level compared to 
HD 19161 2 and HD 108. Its ro tational period is estimated to be 
~64 days ("G runhut et aIJl20lTl) . In the present work, we used a 
set of exposures taken between May 4* and May 9* 2009. 

Finally, we have included t Sco (HD 149438) in our sample. 
Although not an O-type star per se, t Sco is an early B dwarf 
(spectral type B0.2V) just at the limit of the O star mass range. 
Since a magnetic fie ld of th e order 300 G has been clearly de- 
tected on its surface dDonati et al. 2006b), it was justified to add 
it to our otherwise relatively limited sample. The geometry of 
the magnetic field is quite complex, with energy in multipoles 
up to order five. The origin of the field is debated, although a 
fossil origin is preferred. We have used a selection of spectra 
taken between 2005 and 2008 to conduct our analysis. 

To maximize the signal-to-noise ratios, we used the aver- 
age of the observations presented in Table lA.fl for all stars. This 
resulted in spectra with SNR above 1000. 



3. Spectral analysis 

3.1. Model Atmospheres 

The spectroscopic analysis was performed with the atmosphere 
co de CMFGEN (a compl ete description of the code is provided 
bv lHimer(S:Milleiill998l) . The models include stellar winds and 
are computed in spherical geometry. The code solves simulta- 
neously the radiative transfer equation (in the co-moving frame) 
and the rate equations, following an iterative scheme. The com- 
putations provide non-LTE occupation numbers. The tempera- 
ture structure is given by the radiative equilibrium equation and 
the density structure by the mass conservation equation. 

Line-blanketing is included in CMFGEN through a super- 
level formalism in order to reduce the computational cost. Our 
models include H, He, C, N, O, Si, Ne, Mg , Si, S, Fe and Ni. For 
each species, the solar abundances from Grevesse et al. (2007) 
were taken as reference. 

The wind structure is described b y a mass loss rate, a 
terminal velocity and a jS- velocity law (iPauldrach et al.l [19861 
iKudritzki et al.l [l989l) . The subsonic part of the atmosphere is 
described by a p seudo-hydr ostatic structure, computed with the 
TLUSTY code dLanz & Hu benv 2003) with a micro-turbulent 
velocity ^, - 10 km s"'. The velocity structure of the at- 



mosphere is then fully defined by connecting the subsonic 
and the supersonic components. We have used clumped mod- 
els, with a volume filling factor of 0.1 (0.01 for ^ Ori - 
see Bouret et al. (2008.) and HD108 - see lNaze et all (l2008a) ). 
The other wind parameters have been adopted from the litera- 
ture (Donati et al."2006b'; 'B ouret et al.l|2008t iNaze et al.ll2008at 
Grunhut et al. 2009; Wade et al.ll201 idlbl) ! 



Once a model atmosphere is converged, the final synthetic 
spectrum is obtained from a formal solution of the radiative 
transfer equation. A depth-dependent micro-turbulent velocity 
^,(r), equal to 5-10 km s"' at the base of the wind (^o) and to 
about 10 % of the terminal velocity far from the photosphere 
(^oo), is allowed. The value ^, used for the models has very lit- 
tle effect on the atmosphere structure, and consequently on the 
emergent spectrum. On the other hand, may modify the shape 
and strength of the photospheric lines. We found that a value of 
10 km s"' gave satisfactory fits for most stars, except HD 57682 
(the most sharp-lined star of the sample) for which a value of 
5 km s"' was better. 



3.2. Stellar parameters 

The stellar parameters of our best models are gathered in Table|2] 
They correspond to the set of parameters that best fit the ob- 
served spectrum of the respective sample star. The uncertainties 
are indicated for each parameter. We will detail below how we 
determined the uncertainties on the nitrogen abundance. The er- 
rors on effective temperature and surface gravity depend on the 
quality of the fits we obtain, and are used to compute the eiTors 
on luminosity, radius and mass. The main source of uncertainty 
on the luminosity is the distance. 

We derived the effective temperature by the ionization bal- 
ance method. In practice, the relative strength of Hei and 
Hen lines was used. We relied preferentially on He: A4411 
and Hen A4542. Additional diagnostics were: He: /14026, 
He I /14388, Hei ^4712, Hei ^4920, He: i5876. Hen ^4200, 
Hen /15412. In the Of?p stars, the He: lines are partly contami- 
nated by some additional emission the origin of which is unclear. 
The contamination is the largest in the "active" states, i.e. when 
the C III A 4650 lines have their maximum emission strength. 
Consequently, the temperature determination is difficult in those 
states and this explains why we rely on the quiet states whenever 
possible. 

The Balmer line wings are the best surface gravity diagnos- 
tics available in optical spectra. We relied mainly on H/3, Hy and 
H5 to constrain \ogg. In the Of?p stars, the Balmer lines suf- 
fer from the same contamination as the He i lines in the active 
state. Fortunately, the wings, especially the blue one, are rela- 
tively unaffected by this contamination, allowing a safe gravity 
determination. Note that H6 is blended with NIII /14097, render- 
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Figure 1. Example of fit of the nitrogen lines in HD 57682. The observed spectrum is the black solid line. The smooth colored lines 
are models with the best fit parameters and different values of N/H. Thex^ function is evaluated using these models to constrain the 
value of surface nitrogen content. In the present case, N/H = 1.3 x 10""^ is the prefeiTed value (see Fig.|2] red curve). More models 
with additional N/H values have been calculated but only four are shown for clarity. 



ing the use of this line less suitable. It is used as a secondary 
indicator. 

The luminosities are computed in two ways. If calibrated 
UV flux and optical- infrared photometry is available, the spec- 
tral energy distrib ution is fitted. Th e Galactic extinction laws of 
[Siaton (1979) and lHowarthl(ll983h are used to derive the extinc- 
tion. The luminosity is adjusted to match the SED, for a given 
distance. If only optical photometry is available, we derived the 
luminosity from the V band magnit ude, the effective tem pera- 
ture, the bolometric correction (from Martins & Plezll2006) . the 



extinction (calculated from B-V) and the distance. Once known, 
the luminosities directly provide the radii. 

The distance of the Orion stars is set to d = 414 + 50 
pc (Menten et al. 2007). The distance to the Of?p stars is de- 
termined by their sup posed membershi p to OB associations, 
and was taken from iHumphrevsl (Il978h . HD 108 belongs to 
the Cas OB5 assoc iation. We kept the luminosity derived by 
Martins et alJ (120 1 01) assuming a distance of 2.51+0.10 kpc and 
fitting the spectral energy distribution. HD 148937 is a mem- 
ber the Ara OBI association, at a distance of 1.38±0.12 kpc. 
HD 191612 belongs to the Cyg OB3 association, at a distance of 
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2.29±0.12 kpc. T Sco has a reliable Hippacros parallax leading 
to a distance of 145+1 1 pc. 

iGrunhut etaP (l2009h assumed the luminosity of HD 57682 
to be ty pical of a late O-type star and thus adopted the calibrated 
value of Mart ins et alj(l2005h . From this luminosity and the fit to 
the lUE spectrum, they derived a distance of 1 .3 kpc. We adopted 
their luminosity for our own models 

The projected rotational velocities and macroturbulent ve- 
locities have been determined from the fit of He i 4713. A first 
guess of V sin; was obtained from the line Fourier transform. It 
usually led to an upper limit, no zero being observed above the 
noise level. We then used various combinations {V sin/ / Vmac) 
to convolve synthetic spectra and selected the values best repre - 
senting the observed line profile (see also lHowarth et al]l2007h . 
We used a Gaussian profile to take macroturbulent velocities into 
account. 



3.3. Surface nitrogen abundance 

The main purpose of the present work is to derive surface nitro- 
gen abundances as accurately as possible, in order to compare 
them to the predictions from stellar evolution models. 

This task is challenging considering that the optical spec- 
tra of O stars exhibit a small number of nitrogen lines, and that 
most of these lines are relatively weak. In practice, we have used 
all nitrogen lines that our models were able to reproduce. Since 
we have a relatively heteregeneous sample of stars in terms of 
effective temperature, the number of lines used is different for 
each star. The only common lines are the N iii 4510 - 4525 fea- 
tures. These lines are relatively free of contamination and do not 
suffer from severe blending. Their formation do not seem to be 
affected by line-blanketing effects. Their lower levels are con- 
nected to the N III ground level by transitions at 434A, a region 
relatively free of metallic lines. Hence, subtle radiat ive transfer 
effects such as those reported bv lNaiarro et al.l (l2006h are not ex- 
pected. For HD 57682 and t Sco we took advantage of the rela- 
tively large number of N n lines present in the observed spectra. 
For the Of?p stars, we only relied on N iii features since the stars 
are too hot to display N ii lines and N iv 4058 was clearly not 
reproduced. This inconsistency is probably not related to an ion- 
ization issue since for stars with lower Tefj the N iii/N ii balance 
matches the Heii/Hei balances. Rivero Gonzalez et al. (2011b) 
reported that in some of their models, they could not correctly re- 
produce this feature (as well as other N iv lines). Part of the prob- 
lem might be due to t he sensitivity o f these lines to wind prop- 
erties ( Rivero Gonzalez et al.ll201 Ibl) . We did not perform an in- 
vestigation of the wind properties of these objects, and thus we 
refrained from using this line when clear inconsistencies were 
observed. In addition, since magnetic field severely complicates 
the determination of wind properties (due to geometry, clump- 
ing), it is very likely that N iv 4058 is not perfectly predicted by 
our spherical models. Similar c onsiderations can be applied to 
the strong N iii 4634 - 4643 lines. iRivero Gonzalez et al.l(l20lTah 
showed that the onset of the wind was critical in the formation 
of this triplet. We note in addition that blends with O ii lines ex- 
ist at 4643 A. We thus used these lines only when pure absorp- 
tion was observed. Besides, the Of?p stars are known to display 
narrow N m 4634 - 4643 emission lines of unknown origin (e.g. 
iNaze et al.l2008ah . In the case of 0' Ori C, we could reliably use 
only the N m 45 10 - 4525 and N iii 4095 features (the latter being 
located in the wing of H6 and thus subject to more uncertainties). 
All the other N iii features were too weak. 




I — • 1 
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12 3 

N/H [x 10''] 

Figure 2. Chi squares computed for the quiet state of HD 57682 
as a function of the nitrogen abundance. Each color represents 
a parameter set: red for the best fit model (parameters given in 
Tab. |2]i, green for variations of the surface gravity (± 0.1 dex 
with respect to the best fit model) and finally, blue for variations 
of the effective temperature (+ 1000 K with respect to the best fit 
model). The symbols correspond to the models computed. The 
different curves are polynomial fits to the respective set of sym- 
bols (red triangles, blue squares,. . .). The black dot with the error 
bars shows the position of the final N/H value and its uncer- 
tainty. 

Special care was taken to estimate the error bars on the de- 
rived surface nitrogen abundances. To do so, we proceeded as 
follows: for the best fit model (best values of Teff and log^) we 
ran models at different nitrogen abundances. We then computed 
two new series of models with Tgif increased (respectively de- 
creased) by the uncertainty on Teff. We proceeded similarly for 
two series of models with the best estimate of Teff , but log g in- 
creased (decreased) by the uncertainty on logg. We ended up 
with five series of models with 5-8 values of N/H for each se- 
ries. Fig.[T]shows an illustration of the model fits to the observa- 
tions in the case of HD 57682. 

We then estimated the preferred N/H content for each serie. 
For that, we calculated for the entire set of lines 

/=i I 

where M, is the model spectrum, 5,- the observed spectrum, cr,- 
the uncertainty on the observed spectrum and the number of 
wavelength points over which the function is estimated. We de- 
termined the value of N/H that minimizes the chi-square for 
each set of parameters (Teff, \ogg). An illustration of;^'^ versus 
N/H for the various models computed for HD 57682 is given in 
Fig.|2] 

To obtain the final value of N/H, we determined the best fit 
value for each combination (Teff, iogg) and computed the aver- 
age and one cr dispersion of this set of values. The final values 
are given in Table |2] 
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Table 2. Stellar parameters derived from quantitative analysis with atmosphere models. 



Star 


HD191612 


0' OriC 


^ Ori A 


HD57682 


HD108 


HD148937 


r Sco 


Spectral type 


06-8f?p 


07V 


09.7Ib 


09IV 


04-8f?p 


05.5-6f?p 


B0.2V 


Teff (K) 


36000±2000 


38000±2000 


29500+1000 


34500+1000 


35000+2000 


40000+2000 


31000+1000 


logg 


3.75±0.10 


4.20±0.20 


3.25+0.10 


4.00+0.20 


3.50+0.10 


4.00+0.10 


4.10+0.10 


log (L/Lq) 


5.45±0.15 


4.80+0.13 


5.64+0.15 


4.79+0.25 


5.70+0.10 


5.80+0.15 


4.37+0.08 


R. (Rq) 


12.3±1.9 


5.8+1.2 


25.4+1.7 


7.0+0.4 


19.4+2.2 


16.6+1.7 


5.3+0.6 


Mev (Mo) 


37.7Jf 


23.8+3-0 


42 8+'-' 


22.0^/ 


48.8+5-^ 


57 9+^-" 


15 9+° '' 


Mspec (Mq) 


38.7+20.7 


19.7+14.0 


42.0+21.2 


17.8+14.02 


43.3+18.7 


101.0+53.1 


13.1+4.5 


V sin! (km s"') 


< 1 


24 


100 


15 


< 1 


<50 


6 


Vmac (km S"') 


40 


30 


~0 


16 


45 


60 


3 


N/H (number) 


(2.7±1.5) lO-'* 


(5.1+4.5) 10-^ 


(3.3+1.0) 10-' 


(1.3+0.3) 10-^ 


(5.9+4.2) lO-* 


(3.0+1.2) 10-^ 


(1.4+0.2) 10"* 


l2+\og{N/H) 


8.43±0.24 


7.82+0.29 


7.52+0.14 


8.11+0.10 


8.77+0.31 


8.48+0.17 


8.15+0.06 



Notes. The evolutionary masses are determined using the tracks o f iBrottetal .1(1201 Ih with an initial rotational velocity of 300 km s 



To further check the robustness of our determinations, we 
ran a few additional tests. In particular, we tested the effect of a 
change of the microturbulent velocity in the atmosphere model 
computation (not just the output synthetic spectrum). Increasing 
this velocity from 5 to 10 km s"' translates to a decrease of the 
derived abundance by 15-20%. Although not negligible, this is 
within the uncertainties of our determinations. 

Abundance determination have been performed by iMorell 

(1201 ID for two our our sample stars: HD 57682 and r Sco. 
For the latter, our results are in very goo d agreemen t 
(12+log(A^///)=8.15±0.06 versus 8.15±0.20). iPrzvbilla et al.l 
( 1201 Oi) also found a value of 8.16+0.12, in good agreement with 
our finding s. For HP 57 682, we derive a nitrogen content much 
larger than "Morel (|20ll) (8.1 1 ±0.10 versus 7.52±0.25). We at- 
tribute partially this discrepancy to the use of Kurucz models 
in the Morel et al. study. They did not self-consistently derive 
the stellar parameters a nd the abundances, b ut adopted Teff and 
log^ from the study of iGrunhut et al.l ( l2009l) . Another cause of 
the problem can be a different line list used to perform the anal- 
ysis. However, we do not have access to the lines used by Morel 
et al. and cannot test this possibility. 

4. Stellar magnetism and nitrogen enrichment 

4.1. Comparison to predictions of rotating models witiiout 
magnetic field 

In Fig. |3]we show the Hertzsprung-Russel (HR) diagram with 
the positions of the magnetic OB stars analyzed in the present 
st udy. We have cho sen the evolutionary tracks and isochrones 
of iBrott et al.l (1201 ll) . Other evolutionary tracks exist in the lit- 
erature, but the Brott et al. tracks have been computed for a 
large set of initial rotational velocities, allowing an investiga- 
tion of the effects of rotation on the results. In addition, they 
have the same initial metal content as that used in our atmo- 
sphere models, rendering possible comparison between our re- 
sults and the predictions of ev olutionary calculatio ns. Compared 
to the recent Geneva tracks (lEkstrom et al.l I20TTI) they predict 
a stronger nitrogen enrichment at a given age (for moderate and 
zero rotation) during the early evolution representing our sample 
stars. Consequently, our conclusions regarding a possible extra 
enrichment of magnetic stars would be re-enforced if we were 
using the Geneva tracks (see below). In Fig. [3] we have selected 
the evolutionary tracks with an initial rotational velocity of 300 
km s"'. This value is typical since it is supposed to lead to an 
average rotation on the main sequence similar to that observed 
dMaeder & MevnetllloOQ) . 



5.5 



O 



4.5 




4.7 



4.6 4.5 
log (Te„) 



4.4 



Figure 3. Hertzsprung-Russel diagram for the sample magnetic 
OB stars. The iso chrones and the evolutionary tracks are from 
iBrott et al] (1201 ll) and were computed with an initial rotational 
velocity of 300 km s ' . 



From the HR diagram, one can estimate the evolutionary 
mass of each individual star. This is done by simple interpolation 
between the available tracks. The results are given in Table |2] 
The uncertainties reflect the errors on the effective temperature 
and luminosity. The age of the stars can also be tentatively de- 
rived. The ages are just a way of evaluating the degree of evolu- 
tion of a given star They are directly used to build Fig. |4] where 
we show the surface nitrogen content as a function of age. This 
plot, together with the HR diagram, can tell us if a star is evolv- 
ing according to the predictions of standard evolutionary models, 
or if significant deviations are detected. 

Let us take HD 148937 as an example. From the HR dia- 
gram, its evolutionary mass should be 57. Q;^^^'^ Mq. The star is 

^■'5 -0.9 



Myr old. Looking at Fig. 31 we see that for this range of 
ages, the surface nitrogen content of HD 148937 is consistent 
with stars of initial mass 60 Mq. This is thus consistent with the 
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Table 3. Lines used to perform Nitrogen abundance determination. 





HD 191612 


0' Ori C 


4" Ori A 


HD57682 


HD108 


HD148937 


r Sco 


Nil 3995 


- 


- 


V 


V 


- 


- 


V 


Nm4004 


- 


- 


- 


- 


- 


V 


V 


Nn4035 


- 


- 


- 


V 


- 


- 


V 


Nil 4041-4043 


- 


- 


V 


V 


- 


- 


V 


N IV 4058 


- 


- 


- 


V 


- 


- 


V 


N in 4095 


- 


V 


- 


- 


- 


- 


- 


Nin4196 


V 


- 


V 


V 


V 


- 


V 


Nin4216 


V 


- 


- 


V 


V 


- 


V 


Nn4447 


- 


- 


V 


V 


- 


- 


V 


Nm 4510^525 


V 


V 


V 


V 


V 


V 


V 


Nin4537 


V 


- 


- 


- 


V 


V 


V 


Nn4602 


- 


- 


- 


- 


- 


- 


V 


Nil 4607 


- 


- 


- 


V 


- 


- 


V 


Nn4613 


- 


- 


- 


V 


- 


- 


V 


Nn4621 


- 


- 


- 


V 


- 


- 


V 


Nn 4634^643 


- 


- 


- 


V 


- 


- 


V 


Nn4788 








V 






V 


Nn4803 








V 






V 


Nin4907 








V 








Nn4995 








V 






V 


Nn 5000-5011 






V 


V 






V 


Nn5026 














V 


Nn5045 








V 






V 


Nn 5676-5679 






V 


V 






V 



T Sco has an initial mass of about 16 M0. Compared to the 
15 M0 evolutionary tracks, its nitrogen content is rather high. 
Only for an effective temperature at the low limit of the derived 
range (and thus a relatively old age) can the N content be barely 
reconciled with the predictions. HD 57682 also has a nitrogen 
content 2 to 3 times the initial value. This value is marginally 
compatible with a 22 M0 star provided its age is on the upper 
side of the range of possible values. Hence, t Sco and HD 57682 
quite possibly show some extra enrichment. 0'Ori C is very sim- 
ilar in mass but it is too young to have produced a significant 
amount of nitrogen and thus the low A^/// we derive is not unex- 
pected. 

The last star of our sample, ^ Ori A, is the oldest and most 
evolved of the magnetic O stars. At an age of 3.7+0.6 Myr and 
with an initial mass of 42. 8^^ 4, its surface NjH should reach at 
least 2.0 x 10""^. We derive a value about 4 times lower (4.4 + 
1.0 X 10"^). As already noted by Bouret et al. (2008), ^ Ori A is 
thus particularly chemically unevolved for its age. 

So far, we have compared the properties of the magnetic O 
stars to evolutionary models with an initial rotational velocity 
of 300 km s"'. This value leads to an average velocity on the 
main sequence consistent with observations of normal O stars. 
However, one might wonder if this value is appropriate for mag- 
netic stars. The presence of a magnetic field is known to reduce 
the rotation rate as illustrated by Townsend et al. (2010) in the 
case of 0- Ori E. lDonati et al.i (12006a) and Martins et all (1201 Oh 
have also noted the extremely low rotation of HD 191612 and 
HD 108. The question is: were these stars slow rotators already 
on the zero age main sequence? If so, the comparison with low 
y sin; evolutionary tracks is necessary. Table |4]presents the spin 
do wn timescales d ue to magnetic braking calculated according 
to l Ud-Doula et al.l (2009.) . We have also listed the ages esti- 
mate d from Figs. |3] and |5] As noted by iDonati et al.l (l2006ah 
and Mar tins et al.l (12010 ^. HD 191612 and HD 108 have spin 
down timescales marginally compatible with their age, implying 
that they might have experienced a reduction of their rotation 
rate by a factor 2-3. However, given their very low present-day 



0.0004 



z 



0.0002 - 



0.001 : 
o.oooB ; 
60 M„ ° : 

" 0.0004 ; 
0.0002 ; 



2 2.5 3 3.5 4 




Figure 4. A^/// as a function of age. Black solid lines are the pre- 
dictions of evolutionary models including rotational mixing with 
an in itial rotation rate of 300km s ' . Models are from lBrott et al.l 
(1201 II) . The symbols are the magnetic O stars studied in the 
present paper. 



estimated initial mass of HD 148937, from which one can con- 
clude that this star does seem to have normal surface abundance 
for its evolutionary status. A similar conclusion can be drawn for 
HD 191612 and HD 108, although for the latter a hint of extra 
enrichment is not excluded. 
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Table 5. Evolutionary masses derived for different initial rota- 
tional velocities. 



Star 



Mev (300) Mev (100) 



HD 191612 
HD 148937 
HD 108 
0' Ori C 
<rOri A 
HD 57682 

T Sco 



37 1*^-' 

6.2 

57 9+'" 
48.8+^-^ 
23.8+1^ 
42 8+'-'' 

15 9+° '^ 



37 5+^ ' 

+§2 

58.7 12 2 

47.1+'-'^ 

23 6+--^ 

43.0+f' 

21.5+f 

15.7r°-'' 



Notes. The models used are fromHrotteLd] ilOUh . 



V sin/, those stars probably had initial projected rotational veloc- 
ities not higher than a few tens of km s ' . All the other magnetic 
O stars have spin down timescales longer than their age. Given 
their rather low V sin/ (with the exception of ^ Ori A), these 
stars were also most likely slow rotators on the zero age main 
sequence. The youngest star of our sample (0'Ori C) confirms 
this trend: it is located close to the ZAMS and features a low 

V sin/ (by O star standards). 

In view of these comments, we show in Fig. |5] the HR and 
N/H-age diagrams with evolutionary tracks computed for an 
initial rotational velocity of 100 km s"'. Using these plots, we 
can see how the effects of rotation affect the previous results. 
The evolutionary masses derived from these new tracks are given 
in Table H The three Of?p stars (HD 108, HD 148937 and 
HD 191612) all appear to be much more N-rich than predicted 
by the 60 M0 evolutionary track, and thus more enriched than 
expected for their mass, r Sco is now very clearly enriched com- 
pared to the 15 Mq tracks which barely departs from the initial 
N content in the first 6 Myr. For HD 57682, an extra enrichment 
is also clearly present. For 0'Ori C, again the young age pre- 
vents any firm conclusion other than that no dramatic nitrogen 
enrichment is observed. For ( Ori A the surface chemical pat- 
tern is lower than predicted by the evolutionary tracks although 
marginal agreement exists if extreme stellar parameters are used. 

Given the above discussion, the following conclusions can 
be drawn: 



- ^ Ori A: For moderate rotation rate, the low N content cannot 
be explained by standard mixing. Some extra process lead- 
ing to the quenching of rotational mixing is necessary to ex- 
plain the peculiar N/H ratio. The conclusion remains valid 
for low rotation rates, although the agreement with the ex- 
pected N/H ratio is marginal given the uncertainties on the 
age and surface nitrogen content. 

- HD 108, HD 148937, HD 191612, t Sco and HD 57682: 
under the reasonable assumption that these stars were slow 
rotators on the zero main sequence (i.e. initial rotation ve- 
locity of a few tens of km s"'), all stars show a higher de- 
gree of nitrogen enrichment than non-magnetic evolutionary 
models predict. In the less likely case where the initial rota- 
tional velocity was normal (i.e. around 300 km s '), HD 108, 
HD 57682, t Sco and HD 191612 show marginal evidence 
for extra enrichment, while HD 148937 appears normal. 

- ©'Ori C: this star does not show peculiar nitrogen enrich- 
ment whatever the initial rotational speed. This is most likely 
because it is too young. 



4.2. Comparison to non-magnetic stars 

To further investigate the surface chemical properties of the six 
magnetic O stars, we show in Fig. |6] the same HR diagram as in 
Fig. [3] to which we have add ed comparis on O stars. The latter 
objec ts have been taken from lBouret et al.l (2005), Martins et al.l 
(I2OII1) . Bouret etal. (2012, A &A, in prep), Esco lano etal . (2012 , 
A&A, in prep). lPrzvbilla et al. (2010) and Hun ter et al.l ( l2009l) . 
All O stars have been analyzed with CMFGEN in the same way 
as in the present study. In Fig.|6]we encode the value of N/H 
by the size of the symbol for each star (including the magnetic 
O stars). In the upper part of the HR diagram, HD 108 and 
HD 148937 fall in a densely populated group of O supergiants 
with N/H larger than 2.0 x 10""*. They thus appear as similarly 
enriched compared to these objects, and present no peculiarity 
in their chemical pattern provided they were initially rotating at 
about 300 km s"' (but see below). 

HD 191612 is located at a lower luminosity than the two 
other Of?p stars. In that region of the HR diagram, the surface 
nitrogen abundance of the few O stars present a large range of 
values, from less than lO"'* to more than 4.0 x 10 In addition, 
the less evolved stars are more enriched than the more evolved 
objects. Given these puzzling results, the surface abundance of 
HD 191612 do not appear to deviate from the general trend, but 
the large range of N/H values observed in that part of the dia- 
gram prevents any firm conclusion regarding its 'normality'. 

The last star in the upper part of the HRD is ( Ori A. Its 
surface nitrogen content strongly deviates from the values ob- 
served in neighboring stars. The O supergiants analyzed so far 
all show significant N enrichment {N/H > 2.0 x lO""*) while ^ 
Ori A is basically chemically unevolved. This confirms our pre- 
vious conclusions that this star clearly shows peculiar chemical 
properties. 

Going down in the HR diagram, we find O'Ori C and 
HD 57682 very close to each other Only three O stars are lo- 
cated close enough to allow a quantitative comparison. These 
objects show a wide range of N/H values, from less than 1.0 
10 to more than 4.0x 10"^. Both O'Ori C and HD 57682 show 
relatively modest enrichment. This is consistent with the two less 
enriched comparison stars. Hence, the two magnetic stars do not 
seem to stand out as peculiar objects. Finally, close to the loca- 
tion of T Sco we have found very few comparison stars and so 
we cannot draw any relevant conclusion. 

At first sight, this comparison only exacerbates the pecuUar 
properties of ^ Ori A. All the other stars appear similar to non- 
magnetic O stars. There is however an important property of the 
magnetic stars to consider: all (except ( Ori A) are slow rotators. 
As explained above, there are indications that they might have 
been slow rotators already on the ZAMS. As such, they are dif- 
ferent from the comparison stars which have diverse values of 
V sin/ (32 stars with V sin/ from roughly to 330 km s"' with 
an average of 116 km s '). Consequently, finding nitrogen en- 
richment in HD 108 and HD 148937 (the two objects for which 
there are sufficient comparison stars to draw firm conclusions) 
similar to that of comparison stars while they are rotating more 
slowly can be interpretated as the sign of an especially efficient 
mixing inside the star Said differently, the (presumably) initially 
slowly rotating magnetic stars have the same surface N content 
as normal O stars. This is an indirect indication that the presence 
of a magnetic field in those objects tends to favor the transport 
of nitrogen from the core towards the surface. 

Clearly, a larger sample of magnetic and comparison stars 
would be needed to strengthen these conclusions. But the above 
comparisons all indicate that there seems to be two groups of 
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Table 4. Magnetic spin-down timescales and ages. 



Star 


Bp 


M 


R 




M 


''"spin 


age(300) 


age(lOO) 




fkGl 

LiVVJ J 


Mn 


Rn 


km s"^ 


nO"' Mn vr^'l 

LJ.V/ ^'^y J ^ J 


rMvrl 


rMvrl 


rMvrl 


HD 191612 


2.5 


37.7 


12.3 


2400 


0.3 


3.8 


2.8!"- 


2.6+"-*' 
2.0+''-^ 
2 0+^-^ 


HD 148937 


1.0 


57.9 


16.6 


2600 


0.3 


11.3 


2 1 


HD 108 


0.5-2 


48.8 


19.4 


2000 


1.0 


2.0-7.8 


^■^-0.3 


©' Ori C 


1.1 


23.8 


6.5 


2500 


3.0 


3.24 


0+2 ° 

"•"-1.0 


0+2-2 

"■"-1.0 


^Ori A 


0.05-0.1 


42.8 


25.4 


2100 


17.0 


6.5-13.0 


3.6+"'' 


3 6+" '' 


HD 57682 


1.7 


22.0 


7.0 


1200 


0.014 


18.8 


2.3+J-^ 




r Sco 


0.5 


15.9 


5.3 


2000 


0.3 


17.0 


1 3+^-^ 


2 r" 



Notes. The magnetic spin-down timescale is calculated according to lUd-Doula et al.l ( 120091) . Bp is the polar magnetic fiel d. The stellar and win d 
parameters are used to calculate Tspin. We have assumed a value of 0.1 for the moment of inertia constant k (see Eq. 25 of lUd-Doula et al.ll2009h . 
The entry 'age(300)' is the age of the star estimated from the HR diagram built with evolutionary tracks computed for an initial rotational velocity 
of 300 km s"'. Similarly, 'age(lOO)' is the age for V sini = 100 km s"'. 
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Figure 5. Same as Figs. |3] and |4] but for initial rotational velocities of 100 km s ' . The initial masses derived from the evolutionary 
tracks are gathered in Table|5] 



magnetic stars: on one side ^ Ori A which is clearly peculiar, 
with no nitrogen enrichment in spite of its evolved status and 
relatively normal rotational velocity; and the other OB stars, es- 
pecially the Of?p objects, possibly displaying larger enrichment 
than non-magnetic stars provided they were indeed slow rotators 
on the ZAMS. 

4.3. On the nature of magnetic field in O stars 

The origin of magnetic field in massive stars is a matter of de- 
bate. The two main possibilities are either a fossil field or a 
field generated by dynamo processes in the radiative region or 
the convective core. In the former case, the field might come 
from the parental molecular cloud, the convective phases during 
pre-main sequence evolution, of from binarity. Relatively sim- 
ple geometries are expected in that case. For dynamo generated 
magnetic fields, several scenarios exists but the Tayler-Spruit 
dynamo has been the most tested (Spruit 1999, 2002). In this 
case, an initial poloidal field is transfered into a toroidal field 
by differential rotation. Later on, the Pitts-Tayler instability de- 



velops in the toroidal field and gives birth to a poloidal field. 
Consequently, the dynamo loop is closed. Controversy exists at 
present to know if this process is effectively at wor k. In a series 
of ma gneto-hydrodynamical (MHD) simulations, lBraithwa"it3 
(l2006h claim to have observed the onset of the Pitts- Tayler in- 
stability and the dynamo loop. But Z ahn et alj (12007 ') conclude 
differently in their simulations: although the Pitts-Tayler insta- 
bility seems present, it has no effect on the poloidal field which 
slowly decreases due to Ohmic decay. Hence the dynamo loop 
is not closed. A go od summary o f these processes is provided by 
lMaedeii(l2009h and lWalder et'an(l2011h (see also lZahnll201 Ih . 

Can we provide any constraint on the nature of O star mag- 
netism? Mevnet & Maeder (2005) ran evolutionary calculations 
under the assumption that the Tayler-Spruit dynamo was work- 
ing. In that case, they could show that one of the main effects 
was an increase of the mixing and of the efficiency of chemi- 
cal element transport. The reason is the nearly solid body rota- 
tion achieved in the radiative zone in the presence of magnetic 
field. In non magnetic rotating models, chemical elements are 
mainly transported by shear turbulence caused by differential 
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1 , Y< I I I I I I I I I I r 




4.7 4.6 4.5 

log (Te„) 

Figure 6. Same as Fig. [3] with comparis on stars added (gre y 
squares). The comparison stars are from iBouret et al.l (l2005h . 
iMartinsetaT] (1201 ID . Bouret et al. (2012, A&A in prep), 
Escolano et al (2012 , A&A, in prep), iHunter et all (12009 ) and 
iPrzvbilla et al.l (l2010h . The size of the symbols (including the 
six magnetic O stars) is now proportional to the N/H value. The 
evolutionary tracks are for v - 300 km s"' . 

rotation. Meridional circulation - created by thermal imbalance 
due to gravity darkening - contributes little to the transport of 
elements, but is the main way of redistributing angular momen- 
tum. If the star rotates nearly solidly, shear turbulence is reduced, 
but meridional circulation is greatly enhanced so that the net ef- 
fect on the transport of chem ical species is a much stronger mix- 
ing. iMevnet & MaedeJ (12005 ) report an increase of the surface 
N/H of a factor 2 to 5 larger than in non m agnetic rotating mod- 
els for a 15 M© star of age 5 to 10 Mvr. iMaede^ (12009 ) note 
that these results depend on the details of the dynamo process 
at work and on the strength of the magnetic field. However, the 
trend seems to be an significant over-enrichment compared to 
standard rotating models. In t heir more rec e nt inv estigation of 
the role of magnetic braking, Mevnet et alj (1201 lb note that if 
the star is rigidly rotating, braking is so strong that it suppresses 
all surface chemical enrichment. Thus the exact role of magnetic 
field on surface abundances is not clear. 

Our results are meant to provide observational constraints 
that could trigger further theoretical developments. At present, 
they can confront theoretical predictions as follows: 

- dynamo fields: the evolutionary calculations focus on the 
Tayler-Spruit dynamo as explained above. One of the ex- 
pectations of dynamo fields is that, similarly to the sun, they 
show complex morphologies and not simple dipolar struc- 
tures. The constraints on the field geometry of magnetic O 
stars are still scarc e, but ^ Ori A stands out one more time 
dBouret et al.ll2008l) since the first observations indicates the 
possible presence of several magnetic spots. If a dynamo 
process is producing the magnetic field of that star, then 
our results in dicate that no s urface enrichment is observed . 
According to iMevnet et alJ (l201lh - (see also IZahnll201lh 



- this could be the sign of a nearly solid-body rotation and 
strong braking (see above). However, the projected rotational 
velocity of ^ Ori A is not especially low (V sin/ =110 km 
s"'). This is not consistent with a strong braking unless the 
star was initially rotating very fast. But in that case, would 
the high level of rotation be compatible with the absence of 
surface chemical enrichment? Clearly, dedicated simulations 
are necessary to see if a star with a dynamo field can still be 
rotating relatively normally and not have produced any sur- 
face nitrogen excess after several Myrs. 
- fossil fields: it is usually thoug ht that fossil fiel d s hav e 
simple dipolar geometries ( D onati & Landstreefl l2009l) . 
Observations indicate that 0'Ori C, HD 57682, HD 191612 
and HD 148937 have dipolar fields. For HD 108, the phase 
coverage is too poor to provide any constraint, t Sco fea- 
tures a more complex magnetic geometry, but Donati et"an 
(2006b) attributes this to the youth of the star. If a dipole 
should remain after a few Myrs in fossil field magnetic 
stars, higher orders are expected to remain in the early 
evolution, after the formation from the parental molecu- 
lar cloud. Donati et al. (2006b) thus favoured a fossil ori- 
gin for the magnetic field of t Sco. According to Sect. |4] 
all those objects were initially slow rotators on the ZAMS. 
Consequently, the observed surface abundances are higher 
than they would be if there was no magnetic field. In prac- 
tice, the degree of enrichment corresponds to that observed 
in stars with initial V sin/ of 300 km s"', while the initial 
rotational velocity of the magnetic stars was more likely of a 
few tens of km s Unfortunately, there are no simulations 
of the effect of fossil fields on the chemical appearance of 
massive stars to which we can compare these results. We can 
thus simply state that dipolar, probably fossil fields tend to 
lead to extra mixin g in the interio r of O stars. According to 
the results of Mevn etet al.l (1201 ih . this requires the presence 
of strong shear created by differential rotation. Magnetic O 
stars with fossil fields might thus have a steep velocity struc- 
ture. 

Clearly, the effects of magnetic field on the properties of O 
stars are still poorly understood. A larger sample of magnetic 
stars would be needed to see if the apparent separation into two 
groups (slowly rotating enriched stars / normally rotating non- 
enriched stars) is confirmed. Analysis of more comparison stars 
is required too. In parallel, new MHD simulations and evolu- 
tionary models, especially including fossil fields, are necessary 
to better understand the physics at work in the interior and atmo- 
sphere of magnetic O stars. 

5. Summary and conclusions 

We have conducted a quantitative spectroscopic analysis of the 
six magnetic O stars known to date together with the early B 
star T Sco. Spectra collected with ESPADONS at CFHT and 
NARVAL at TBL have been used. Atmosphere models com- 
puted with the code CMFGEN have been computed. We have 
determined the stellar parameters and the surface nitrogen abun- 
dance. A careful investigation of the uncertainties associated 
with N/H has been pursued. 

Our results can be divided into two parts: 

• The surface N/H of ( Ori A - the star with the weakest 
magnetic field - is lower than the predictions of evolution- 
ary models. For that star the presence of the magnetic field 
probably quenches the transport of chemical elements. The 
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conditions under which such a process can happen are not 
clear. Evolutionary models accounting for magnetic braking 
predict no chemical enrichment in solid body rotating stars. 
But their rotational velocities are also rather low, whereas 
C, Ori A appears to have a rather standard V sin;. 
• The other magnetic O stars, especially the three Of?p stars, 
display surface nitrogen content consistent with those of 
non-magnetic O stars with initial V sin? of 300 km s"'. 
However, these magnetic O stars were most likely slow rota- 
tor on the zero age main sequence (with values of Y sin; of a 
few tens of km s"' at most). Consequently, they feature larger 
surface iV/// than they should if they did not host magnetic 
fields as also indicated by their nitrogen excess compared to 
rotating evolutionary models. 

The presence of a magnetic field in O stars thus appears 
to produce two types of chemical peculiarities: either an extra 
enrichment, or the absence of enrichment. Confrontations with 
evolutionary predictions in the case of dynamo fields indicate 
that solid body rotation might explain the properties of the non- 
enriched star, but its relatively normal rotational velocity is a 
puzzle in that context. For the N-rich magnetic O stars, a fos- 
sil field is usually assumed. Our results should be a basis for 
comparison with future evolutionary calculations including fos- 
sil fields since such evolutionary models do not exist at present. 
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Appendix A: Journal of observations 

Table lA. ll provides the complete journal of observations. 
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Table A.l. Journal of observations. 



Star 


Inslrumcnt 


date 


UT 


Exp. lime (sec) 


S/N 


HD 108 


NARVAL 


Jul 2P' 2009 


00:52:37 


4 X 750 


150-350 




NARVAL 


Jul 21" 2009 


02:18:22 


4x750 


200-500 




NARVAL 


Jul 2P' 2009 


03:11:49 


4 X 750 


200-500 




NARVAL 


Jul 25* 2009 


00:41:56 


4 X 750 


200-500 




NARVAL 


Jul 25* 2009 


01:35:23 


4 X 750 


200-500 




NARVAL 


Jul 25* 2009 


02:28:49 


4 X 750 


200-500 




NARVAL 


Jul 26" 2009 


01:47:60 


4 X 700 


100-350 




NARVAL 


Jul 26" 2009 


02:38:05 


4 X 700 


100-400 




NARVAL 


Jul 26" 2009 


03:28:11 


4 X 700 


150-450 




NARVAL 


Jul 27* 2009 


00:53:23 


4 X 806 


100-300 




NARVAL 


Jul 27* 2009 


01:50:33 


4 X 806 


120-370 




NARVAL 


Jul 27* 2009 


02:47:43 


4 X 806 


100-350 




NARVAL 


Jul 28* 2009 


00:51:16 


4 X 900 


100-300 




NARVAL 


Jul 28" 2009 


01:54:42 


4 X 900 


100-400 




NARVAL 


Jul 28* 2009 


02:58:06 


4 X 900 


150-450 




NARVAL 


Jul 29* 2009 


00:39:56 


4x 900 


200-550 




NARVAL 


Jul 29* 2009 


01:43:21 


4x 900 


250-600 




NARVAL 


Jul 29* 2009 


02:46:46 


4x 900 


200-500 




NARVAL 


Jul 30* 2009 


00:56:23 


4x 900 


250-630 




NARVAL 


Jul 30* 2009 


01:59:50 


4x 900 


250-600 




NARVAL 


Jul 30* 2009 


03:03:16 


4x900 


250-600 




NARVAL 


Jul 31" 2009 


00:30:28 


4x812 


200-500 




NARVAL 


Jul 31" 2009 


01:28:02 


4x 812 


200-500 




NARVAL 


Jul3P' 2009 


02:25:37 


4x 812 


200-550 




NARVAL 


Jul3P' 2009 


03:23:11 


4x 812 


200-500 




NARVAL 


Aug P' 2009 


00:55:57 


4 x 812 


200-550 




NARVAL 


Aug 1" 2009 


01:57:06 


4 X 812 


250-600 




NARVAL 


Aug P' 2009 


02:54:39 


4 X 812 


200-470 




NARVAL 


Aug 2"^ 2009 


00:50:14 


4 X 900 


200-470 




JNAKVAL, 


Aug 2 /U09 




4 X you 






NARVAL 


Aug 2'"' 2009 


02:57:24 


4 X 900 


150-400 




NARVAL 


Aug 4" 2009 


01:24:42 


4 X 900 


200-470 




NARVAL 


Aug 4* 2009 


02:28:07 


4 X 900 


250-600 




NARVAL 


Aug 4* 2009 


03:31:34 


4 X 900 


200-570 




NARVAL 


Aug 5* 2009 


00:22:43 


4 X 900 


200-500 




M A D\?A T 

iNAKVAL 


Aug J zuuy 


ui.zo.uy 


4 X yuu 


zUU— OUU 




TVT A r*T TAT 

NARVAL 


Aug 5" 2009 


02:29:35 


4 X 900 


200-600 


HD 148937 


ESR\DONS 


Jun 20" 2010 


09:10:53 


4 X 600 


300-850 




ESPADONS 


Jun 21 2010 


AO O 1 

08:31:20 


4 X 600 


250-750 




ESPADONS 


Jun 22 2010 


09: 9:60 


4x 600 


170-550 




ESPADONS 


T ril 1 r\ 

Jun 23 2010 


08:11:34 


4 X 600 


1 C A A f\f\ 

150-400 




ESPADONS 


T Ath '^nt n 

Jun 24 2010 


08: 8:21 


4 X 600 


O AA OAA 

300-800 




ESPADONS 


Jun 25" 2010 


08:16:39 


4 X 600 


C A O A A 

3250-800 






Anrr 011^ OHrtQ 

Aug ZUUo 




A v inn 

4 X jUU 


iUlf— jUU 




ESPADONS 


Aug 23"' 2008 


10:04:20 


4x 300 


100-350 




ESPADONS 


Aug 23'^'' 2008 


10:25:05 


4x 300 


100-350 




ESPADONS 


Aug 23'^'' 2008 


10:45:49 


4x 300 


100-350 




ESPADONS 


Aug 23'^'' 2008 


11:07:00 


4x300 


100-350 




ESPADONS 


Aug 23'''' 2008 


11:27:45 


4x300 


100-350 




ESPADONS 


Aug 23"* 2008 


11:48:29 


4x300 


100-350 




ESPADONS 


Aug 23'-'' 2008 


12:09:13 


4x300 


100-350 
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Table A.l. Continued. 



Star 


Tnctmmi^Tif 
Ills LI LllllCllL 




ITT 


XjAU. LllllC l^SCL^ 


S/N 

Oj IN 


(5)1 OH C 


NARVAT 


Mar in"' 7nn7 




A V 400 






NARVAL 


Mar 10* 2007 


08-5901 


4 X 400 


450-1250 


/ rVri A 




Cini 1 9007 


00- ^^9-^^ 


AQ V A V 90 


780 QQO 




MAR VAT 


Cini 1 Qth 9007 




Q V ZL V ZlO 
O A A H-U 


1010 10S0 




TSJARVAT 


Cini 90th 9007 
WLL ZU ZUU / 


yjyj.jy.jy 


A A V Zl V ZlO 


1990 1470 




MAR VAT 


Orf 9 1 ^t 9007 

WLL A i 1 


9^-OfS- 1 S 

Z.J>.V7U. i J 


ZLS V ZL V ZLO 

H-O A H- A H-U 


810 1 4(S0 

iU — iH-UU 




NARVAL 


Oct 22"'' 2007 


23:45:52 


48 X 4 X 40 


1090 ~ 1480 




NARVAI 


Oct IV^ 7007 


23:13:45 


AS V 4 Y 40 


1 030 - 1 480 




NARVAT 


Ort 24* 2007 




49 Y 4 Y 40 


1 200 - 1 470 

lZ,Uv/ — iT-ZU 


nu J / uoz. 




Mav 9000 




4 V (SOO 

H- A UUU 


SOO 1 1 ^0 

J UU — i 1 J u 






Mav 900Q 


0(S-99-17 


4 X 540 


4^0 1 000 

T-JU — iUUU 






\/\ct\T 7th 900Q 


UD.H-D.H- i 


4 V '^40 


■^SO 700 
J JU— /UU 




ESPADONS 


May 8* 2009 


06:21:09 


4x 540 


350-800 




ESPADONS 


May 9* 2009 


06:20:36 


4x540 


300-600 


T ScO 


ESPADONS 


May 23'" 2005 


09:14:25 


4x60 


700-1550 




ESPADONS 


Sep 20* 2005 


05:01:56 


4x30 


600-1200 




ESPADONS 


Aug 5* 2006 


05:25:22 


4x 30 


500-1150 




ESPADONS 


Mar 4* 2007 


16:23:15 


4x30 


700-1350 




ESPADONS 


Jun 29* 2008 


08:03:07 


4x45 


400-1050 



